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THE CYLINDER AMD SEMICYLINDER IN SUBSONIC, FLOW 


riNTRODUCTION 

In studying the diffraction of shock waves around various two- 
dimensional obstacles we have observed that flow separation and the for- 
mation of vortices contributes in an important way to transient loading 
of the obstacle. The cases of a cylinder and semicylinder are especially 
interesting because the breakaway point is not clearly defined as it is 
for objects having sharp corners. Accordingly a nvunber of experiments 
have been made in the shock tube to observe the influence of Reynolds 
number and Mach number on the transient flow patterns about a cylinder and 
about a semicylinder mounted on a smooth plane. Some differences might be 
anticipated since the pleuae would impose a symmetry on the flow and pro- 
duce a viscous boundary layer for which there is no counterpart with the 
cylinder. 


In the course of these experiments it was noted that a condition 
of steady subsonic flow about both the cylinder and semicylinder was ap- 
proached. Thus a comparison with von Kannan's^»^ theoretical calculation 
of the drag on a cylinder, from certain characteristics of its wake or 
"vortex street", was underteiken. 


EXPERIMENTAL ARRANGEMENT 

The theory and operation of the shock tube along with the tech- 
nique of obtaining the density and pressure flel^ arotind a model have been 
well described by BleaJmey, Weimer and Fletcher.-^ In brief a shock tube 
consists of regions of high and low pressure separated by a cellophane 
partition. ■ When this diaphragm is pionctured, a compression wave which 
rapidly steepens into a shock travels down the lower pressure section of 
the tube (the channel) and a rarefaction propagates up the higher pressure 
section (the chamber). The density field aroiand the model due to the pas- 
sage of the shock wave and the air flowing along behind it is determined 
by use of a Mach-Zehnder interferometer. For two dimensional flow the den- 
sity change is simply proportional to the fringe shift. 

The Mach number of the shock wave depends only on the initial 
pressure ratio across the diaphreigm (if the initial temperatures are the 
same throughout the tube) . The Mach number, pressure and density of the 
flow behind the shock wave (the conditions "infinitely" far from the 
model) are determined by the initial conditions of the gas throiagh which 
the shock travels and the Mach niimber of the shock. The Mach number of 
the flow at a given point in the test section remains constant until the 
incident shock returns after reflection from the end of the tube. 

The models used for most of the experiments reported here were 
a half- inch diameter polished brass cylinder squeezed between viewing 
windows by a rubber gasket and an equal diameter semicylinder mounted 
on a flat steel plate extending from the rear wall of the test section 
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about two feet behind the viewing windows to a point one foot upstreaj^. 

The Reynolds number Re^ based on cylinder diameter, varied between 10 and 
10^ for the range of flow speeds investigated. The time between the 
arrival of the incident shock and its return from the end of the tube was 
l400 microseconds for the fastest flow, MO. 735 and 2200 microseconds for 
the lowest Mach number, O.I56. A dimensionless unit of time X. is defined 
as the time elapsed since the shock first struck the cylindgjr.jcLlvlded • by 
the -time to travel one cylinder diameter. Some tests were made with a 
2 1/2" diameter cylinder but accurate results were limited to early flow 
times by the top and bottom wall interaction. For the I/2". models no. 
such interaction was detected. Plate I is a typical interferogram show- 
ing the disturbed fringe pattern just after the Incident, shock lias passed 
the cylinder. 


EQUATIONS AND SYMBOLS 


To determine the density p at any point in the field from the 


measured fringe shift between t^e initial condition A and final 
condition B and the initial density .. we use 


= 1 + O.OI 

where is the density at standard temperature and pressure and COI86 

is a constant determined from the rate at which fringes pass a given 
point in the field of view as the pressure in the tube is varied,' 

The pressure p is related td the density by 



where f" is the ratio of specific heats, (for air = 7/5 =1?^) and the 
subscript Ao refers to the stagnation pressure corresponding to the stAte A, 

From the pressure field we pick out for particular attention tho 
quantities; 

W _ ^3 ”^1 tJie overpressure on the cylinder in units 

Pi 

of the initial pressure in the test section p , Subscript 1 denotes 
initial conditions in the test section, the subscript 2 refers to the 
conditions just behind the undisturbed incident shock, and subscript 3 
refers to the conditions at a point in the field after ^the passage of 
the shock. 
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the pressure coefficient. 




is the Mach number of the flow behind the incident shock. 



© d © 


the drag coefficient per 
unit length based on 
cylinder diameter. 


For each of several sets of initial conditions we determined 
the variation of these quantities with X • Figures la to 4f are repre- 
sentative of the density fields at various times for the initial condi- 
tions given on the figures. In all cases the Incident shock travels to 
the leftt The contours are lines of constant fringe shift 8si . is 

the incident shocks the shock reflected from the cylinder and ^Sg 

the portion of the incident shock which sweeps back around the cylinder. 

SSo< is the slipstream growing from the intersection of S_ and , 

X 'i\ 

the "Mach stem" or altered incident shock, SS ^similarly grows from 
the triple point of S„ , S„ (the "secondary Mach stem") and S„, 

M- iYl„ p 


* By accident the convention adopted in this report is just the opposite 
of that used in our other technical papers. 


EXPERIMENTAL RESULTS 

1. Cylinder with ^2 _ 1.25 , M = .156 , Re =32,200 

■n • ^ ■ 


Representative density fields around the cylinder for values 
of 't from 2.3 to 37*2 for a shock of pressure ratio 1.25 are given in 
Figures la - le. The variation of Cp (and W) with angle for these same 
times is shown in Graphs la 7 le respectively, In addition curves are 
included which exhibit the time dependence of various quantities 
measured from the series of Figures and Graphs. In these figures and 
graphs one can easily follow the decreasing pressure on the front 
(right) side of the cylinder from the initial head on reflection value 
as the reflected wave (Sp ) goes out, till by Figure le a rather close 

approximation to steady potential flow is in evidence. The initial 
high pressure on the rear due to passage of the two branches of 
the original shock decays more slowly than on the front so that 
the drag is for a short time negative (Figure lb) . Note that at 
this time the field is much like potential flow except that the 
pressures are uniformly higher. In. Figure Ic we have dips in the 


3 


pressure at points A. and B .showing the beginnings of separation of flow and 
the formation of vortices. In the later pictures the growth and rearward 
movement of these symmetrically placed vortices can be clearly traced as 
they reduce the earlier high rear pressure and increase the drag . In the 
latest picture they have moved far enough rearward for the pressure to climb 
again and reduce the drag. At this stage the flow separates at 88 and the 
vortices are still moving back symmetrically. We expect that were we able 
to follow them longer we would discover a regular' alternation of vortices 
similar to the results at ^2 — 1«8 below. 

>1 

2o Cylinder at ^2 1.8 , M Re -=76,600 


Figures and Graphs 2a ^ 2m give our data at this shock strength. 
One can note again in this series a decay of the front pressure from' the 
head on reflection value to approximately the value for potential flow. 
Comparisons of theise results with those of the previous section show that 
at the higher strength^ the shock sweeps more slowly around toward the 
front of the cylinder, and the vortices form more rapidly, in fact right 
behind the sweeping shock. The two sets of slipstreams SS*^ and SS^ are 
also clearly, visible . Note that the more rapid formation of vortices 
causes the pressure on the rear to d,ecay much more rapidly than, on the 
front so that on Curve 2 .the drag never becpmes negative. In fact, as may 
be seen from Curve 2 the drag reaches a maximum between Figures 2c and 2d 
where the vortices are growing large but are still right on the rear of the 
cylinder. As they mpve rearward (through Figure 2g) the drag drops since 
the pressure on the rear of the cylinder rises and is still decreasing on 
the front. Note in Graph 2c the small inflections at points A and B, about 
100^. Perhaps these points are where the flow begins to separate from the 
cylinder surface* In Figure 2e, note the high pressure regioct where the 
streamlines from the two vortices meet head on. This persists through 
Figure 2g ^and similar high pressure areas may be seen later among the 
successive vortices. By Figure 2f Jhe separation appears to have moved 
forward to points A and B, about 82 from the nose, which is consistent 
with the fact that the vortices have grown larger (though as they mpve 
rearward they become less deep, their depth decaying approximately ex- 
ponentially with the distance from the cylinder center) , In Figure 2g 
we can see a phenomenon unobserved in the time we had to follow the 
vortices produced by a 1.25 shock namely one vortex gets ^head and 
thus the fplrce on the cylinder becomes unsymmetrical (the closer vortex 
exerting the greater influence). We note that ih spite of the less enihg 
influence of this farther vortex that the drag is again on the increase. 
This is probably due both to the beginning of turbulence preceding the 
formation of a third vortex, and to the continued growth of the second 
vortex (close one to the cylinder) . . Also its movement closer to the 
axis, causes it to exert all its influence in the horizontal direction 
rather than obliquely. 

^evenX later in Figure 2b we have a third vortex well formed 
and growing, the drag increasing still further because of it. The 
separation point of the flow is not apparent pn the top, but on the 
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bottom it has remained at about 83°. From here on (Figures 2i - 2m) the 
vortices are shed alternately from the top and bottom of the cylinder form- 
ing examples of the Karman vortex street. As the vortices form and leave, 
the point of separation seems to snap back and forth, though its movement 
is difficult to follov clearly. The average pressure drag coefficient re- 
mains at about I.5 from 'T7 - 30 

An example of the origin of fluctuations observed at later 
times is shown in Figure 2k. The fifth vortex has moved down almost to 
the axis so that the rear pressure has dropped far below its average 
value, causing a sharp peak in the drag curve. All the data available 
are consistent with the assumption that similar maxima in drag are asso-r 
elated with each vortex leaving the cylinder so that it is subject to a 
periodic horizontal and vertical force. Curve 2a has been drawn on this 
basis even though insufficient points are available to fully Justify such 
a procedure. 

The pressure at the nose and at 45° intervals around the 
cylinder is plotted as a function of time in Curves 2o^j^, and 2o^^ with a 
detail of early tlines in Curve 20 /^ initial values are calculated 

as follows; for the front from the pressure due to reflection of the 
shock; for 45 , 90°, 135 , using the appropriate pictures in Reference 
4 (No. 508, No. 405 Vj for 180° from measurement of S on a picture taken 
especially for the purpose. The pressure on the front may be seen to 
approach approximately the stagnation pressure value. At 45° the pressure 
is about 30?^ too high for potential flow, at 90° about two and a half 
times as high as potential flow, at 135° about 155^ lower • and l80° about 
half the potential flow value. The discontinuities in Curves 2«^, - 2e^j 
(e.g., at T 2,2 on the 135° Curve) are due to the passage of the sweepr 
ing shock Sg. Its decay as it moves forward from the rear of the cylinder 
may be followed in the several curves. The fluctuations in pressure on 
the cylinder rear due to the formation and shedding of vortices may be 
clearly followed in the curves for 135° and l80° and to some extent for 
90 . We will consider these curves further in comparison with the semi- 
cylinder results. 

Similar results with respect to drag, pressure curves and vor- 
tex street were obtained for some series shot at approximately the same 
Mach number but different Reynolds numbers. A few of these results are 
marked (+) on Curve 2a, 

The existence of wjl^ formed vortex streets suggested compari- 
son with von Karman* s theory ^ . An analysis of the stability of various 
conceivable configurations of vortex trails led von Karman to the conclu- 
sion that the only stable configuration is the assymmetrical type with 


He also found an expression for the drag D on the cylinder per unit 
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length in terms of the strength of the vortices P » thfe density P , the 
free stream velocity the vortex velocity relative to the stream u. ^ 

r' 

^ flTi 

r 


A 


we can reduce the above expression to 



Measurement of Figures 2a - 2n and others not included here 
gives l^/XsO .310 for all yortices and jk/Jt - 0.272 for those heyond about 
one diameter from the cylinder. This is in good agreement with the 
theory. The ratio of vortex separation to cylinder diameter is 3.1^* 

To find a drag coefficient from the theory an estimate must he 
made of the speed of the vortex centers as well. The points on Curve 2^ 
represent the positions of the various vortices appearing in each picture 
plotted against the <x of that. picture. For example, the three vortices 
of Figure 2h are represented by points above T= 28.65 at the height cor-r 
responding to thsir respective distances from the cylinder center. Since 
each picture represents a separate experiment in the shock tube it is 
sometimes difficult to tell whether the tpp or bottom vortex left the 
cylinder first. All of the Figures 2a - 2m have been oriented so that 
the first vortex left the top according to our best guess. The numbers 
on the drag maxima in Curve 2a as measured from the. pressure are assigned 
accordingly. 

The path of each vortex may be traced fairly unambiguously in 
Curve 2 fi . . An idea of the reproducibility of the experiments may be 

seen in bhe two sets of points at *X=-62 which come from separate firings 
of the tube. The four symbols O p 6 represent series taken under 
slightly , different conditions of temperature and humidity, 

Values for the vortex velocity can be obtained from Curve 2 
The average slope of the lines drawn gives u.^ — 0.2to times the shock 
velocity. Since the velocity V of the fre.e^ stream with respect to the 
body is 0,338 times the shock velocity we f ind _ V _;290. Sub- 
stitution in the formula gives Cn ■= 1.28 whichys in fair agreement 
with the value 1.5 determined from an average of Curve 2a. 

A periodic lift force also acts on the cylinder as vortices 
are shed. This force has half the frequency of the fluctuations in 
drag since every other vortex leaves a given side of the cylinder. 
Sufficient experimental points have not been obtained to trace through 


and the ratio Jk/ i, 


Where 


Using 


V = f r^(v - 

. . XL ^ 

^ ~ XI i 


c, = 


J.O j -\r ‘ 
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this variation completely hut all those available are consistent with the 
assignment of vortices shown in Curve 2a. Values of lift coefficient 
vs. T are plotted in Curve 2b, The solid line should not be taken 
seriously, it is only drawn as a qualitative suggestion. During early 
stages the flow is symmetrical and the lift negligible. After T = 20 
oscillatipns build up to magnitudes of 0,9« So far as we know this 
has not been measured before because the shedding frequency is too fast 
for response by a wind tunnel lift balance. Measurement with a strain 
gauge would presumably yield results on the fluctuations as well as on the 
average values . 


3. Semicylinder at ^ =1,8, M =.0,4, Re -=76,000. 

Pi 


As might have been expected, the density field around a semlr 
cylinder is essentially the same as for a cylinder during the early stages 
of flow. Even when the flow about the cylinder has become unsymmetrical, 
(Figures 3d, 3©) certain portions of the patterns are qualitatively similar, 
including the high pressure region A where a stagnation point, forms 
on the plate in somewhat the way one forms behind the two vortices in Figures 
2e, 2f , The corner between the plate and the semicylinder rear forms a 
better trap for still air than the unprotected rear of the cylinder, 
however, so the pressure remains systematically higher behini the semi- 
cylinder than behind the cylinder. Thus the measured drag on the semi-; 
cylinder shown in Curve 3 lias an average value C»=.0 -9o compared to Ct = 1.5 
for the cylinder. 

At the later times the pattern is no longer always even quali- 
tatively similar. For example, in Figure 3f^ T = 23-4 the second vortex is 
already well formed while the first is apparently beginning to break up. 
Figure 2g for the cylinder, ( T=21,5) > on the other hand shows both 
vortices still well formed moving together, and little sign of a third 
vortex (which would correspond to the second on the semicylinder) . In 
Figure 3S ( X-30.2) the second vortex has grown larger and moved back 
allowing the dreig to further decrease. The positions of the first and 
second vortices in this picture correspond roughly to the first and 
second from the upper surface of the cylinder (i.e., the first and third) 
in Figure 2h ( T=28.65). The pressure distributions over the rear in 
Graphs 3S and 2h are quite different because the vortices cannot grow as 
large behind the semicylinder before being swept away as they can behind 
the full cylinder. 

In Figure 3 I 1 we have what seems to be the third vortex still 
growing and then another meiximum in the C-o curve. By Figure 31 the third 
vortex has pretty well gone so the drag is again low, In these pictures 
the position of the point of separation of the flow from the cylinder 
‘ surface seems to be furtlaer back than in set 2, in accord with the smalleV’ 
vortices. No separation point further forward than 86 was noted. 

That three vortices form by T = 50 would be suspected from 
the behavior of the whole cylinder, from which the sixth vortex is de- 
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veloping by T = 50- That each vortex would form On the semicylinder some-! 
what sooner than the corresponding vortex on, say, the top of the whole 
cylinder is also to be expected since there would be no vortex from the 
bottom to interfere. Understandable also is the more rapid breakup of the 
semicylinder vortices due to the influence of the plate. 

Some insight into the origin of the differences in drag be- 
tween the cylinder and semicylinder may be gained by studying the varia- 
tion of pressure with time for the locations shown in Curves 2«<'i - 
The pressure on the front is seen to approach the stagnation pressure 
about as fast in both models even though a boundary layer forms along the 
plate ahead of the semicylinder. At 45 the pressure is uniformly the 
same in the t^o models over the time range studied within experimental 
error. At 90° a slightly higher pressure appears after x -20 or so. 

This probably arises because of the smaller size of the vortices shed 
from the semicylinder. At 135°^ though again the differences before 
•Z -20 are not significant, the later times show clearly that the pres- 
sure on the semicylinder is always higher than on the whole cylinder. 

This again is because the vortices are not so large, i.e., because there 
is a region of dead air which the plate protects from being swept away. 

The curve for l80° shows clearest of all the presence of this high pres- 
sure region. There seems also to be a difference between the pressure 
on the rears of the two models between T-5 and 10, possibly because of 
the difference in the manner of growth of the first vortices in the two 
cases. 

4. Cylinder at ^_3,0 , M =.735, Re ■= 30,800 


As the flow Mach number is increased, more and more phenomena 
resulting from the compressibility of the air begin to appear. Some 
highly Interesting effects adjacent to regions of locally supersonic 
flow may be seen in Figures 4a - 4d. 

In Figure 4a for example we note that the sweeping shock seems 
to be having difficulty going back upstream around the cylinder, having 
reached only some 107° in x "=■ 4.76 (compare Figure 2c where in X “ 5°1 
the shock is just about at the cylinder nose, some 95° further) . The 
vortices, too are forming more rapidly and have developed little shock 
waves like horns on their sides away from the axis. The slip streams are 
again visible along with the reflected and sweeping shocks. 

The next picture (Figure 4b) shows the sweep shock having given 
up trying, in fact to have been partly swept away itself (S^ between the 
two vortices) and the rest of it curled past the cylinder. The drag has 
dropped nearly 405& due to the decay of the reflection high in the front 
and the rapid removal of the first two vortices (which still have their . 
horns) . Figure 4c shows the remnants of the sweep shock still being taken 
slowly back by the flow, while two new shocks (A and B on graph and fig- 
ure) appear about t 75° from the cylinder nose. The vortices with their 
horns still move back, the rear pressure continues to rise and the drag 


8 


to drop, 


Later, we find the sweep shock disappears entirely and regular 
alternation of vortices takes, place , These vortices have horns till they 
are a diameter or two hack of the cylinder, hy which time they have 
widened put so. that they no longer have regions of supersonic rotational 
velocity and the horns have disappeared. An example of the later appear- 
ance of the flow is. Figure Ud which shows the assymmetrical character of 
this flow due to the formation of vortices. We did not take enough 
pictures to follow the drag fluctuations in detail or to check the Karman 
Ci> prediction. The points we have form a smooth curve suggesting that the 
fluctuations in the dreig with time are not large, The dreig coefficient, 
as may be seen from Curve k, rises as the later vortices are formed and 
shed and reaches a value of 2.2 at x — 60. From what pictures we have 
we found that iv/JL .276 (including all vortices) suggesting, perhaps that 
with increasing Mach number decreases. 

In all of the analysis a tacit assumption has been made that 
the flow patterns observed are accurately two-dlmentional. Previous 
experience supports this idea since the boundary layer which forms on the 
side walls of the tube is very thin compared to the 4” tube width, Be- 
cause the shock and slipstreams observed in this group of experiments 
appear sharp in the pictures it also seems reasonable to assume that the 
flow past the 1/2" diameter cylinder is actually two-dimensional over the 
4" span. 


CONCnJSIONS 

Our expectation that the flat plate would affect the transient 
and steady flow pattern around a cylindrical obstacle was borne out by 
experiments with shock waves of^^/^,=1.8. In particular it was found 

that the rearward extending plate protected a high pressure "dead air" 
region behind the cylinder, which lowered considerably its drag as com- 
pared to a free cylinder. 

Fluctuations in the pressure drag coefficient with the forma- 
tion and shedding of vortices alternately from the top and bottom of the 
free cylinder were found to be on the order of - 205^ from the meeoi ralue 
of 1,5 (for Mach number 0.40, Reynolds number 76,000) . Similar results 
were also obtained for Reynolds numbers from 12,000 to 86,000 at approxi- 
mately this Mach number. Our data suggest but are insufficient to prove, 
that there is little variation of the average drag coefficient within this 
range of Reynolds numbers for Mach number 0.40. The lift coefficient 
varies over about - 0.9 with half the frequency of the drag fluctuations. 

The presence of the plate upstream from the semicylinder was 
found to have little detectable effect. Downstream, however, somewhat 
higher pressures were observed and the drag coefficient was reduced to 
about 0,90. 

von Karman' s theory for ithe stable configuration of the wake 
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of vortices behind, the cylinder was approximately verified, not only for 
Mach numbers around OA, Reynolds numbers in the range 10v«» 10^ but for a 
few cases of Mach number as high as 0,73 (Reynolds number 31>000) . 

At M ,l|-0 JL/H was found to be .272 as compared with a predicted value 

. 281 . 

Calculation of the drag coefficient from the character of the 
vortex street according to von Karman’ s theory was found to give a value 
for Cp “a 1.28. No wind tunnel experiments of which we are aware have been 
done for these conditions. Several investigators have studied this 
range of Reynolds number at very low Mach numbers and eigree onCp-l.?, 

Supporting experiments at Mach numbers 0.l6 revealed flow pat^ 
terns qualitatively similar to potential flow, but with the pressure 
uniformly higher. These flow patterns persist for only a few shock cross- 
irig times before the flow begins to separate and form vortices. At higher 
Mach numbers the vortices form so quickly that a pattern similar to po- 
tbntial theory is never observed. The vortex center depth was found to 
decay approximately exponentially with distance from the cylinder center, 

The experiments at Mach number 0.73 revealed complex interac- 
tions between the flow and the various diffracted shock waves, The pres- 
ence of locally supersonic flow was indicated by shock waves attached to 
the cylinder surfaces and to vortices close to the cylinder, 
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